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ABSTRACT 


The measurement of the human describing functior and. 
remnant in a compensatory tracking task. is undertaken.. 
These measurements are obtained through the application. of. 
the fast Fourier transform technique on a hybrid (analog— 
digital) computer. This method processes the data in. real 
time with minimal core storage and the results are available: 


immediately upon completion of the tracking run.. 
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I. INTRODUCTION 


А. BACKGROUND 

Many of the important tasks performed by pilots are- akin 
to those performed by linear servomechanisms. In. situations- 
such as this, the pilot can be modeled by a set. of constant-- 
coefficient linear differential equations... In the frequency 
domain, such a model is often referred to as a “human. pilot: 
describing function." The term "describing function" is: 
preferred to "transfer function" to emphasize the fact. that 
this pilot model is approximating a nonlinear element. and: 
is valid only for the particular inputs, system dynamics: and 
task at hand. 

The pilot-describing function is useful in studying two 
classes of problems. First, the describing functions 
measured in the piloted simulation of a. given aircraft. 
and task can be utilized in the subsequent stability and 
control analysis of this aircraft. Once the pilot's 
describing function for a particular task has been. measured, 
he can be analytically replaced by his describing function 
in the analyses normally associated with the study of linear 
feedback systems. Second, actual measurement of pilot 
describing functions in ground simulation or: in. flight. tests: 
can be used to determine how a particular aircraft or flight. 
task affects pilot behavior. Knowledge of pilot describing 


functions consequently provides valuable information for 





the aircraft designer. 

It is the problem of human describing function.measure- 
ment which forms the basis of this research. The hybrid 
computer program which has resulted can be utilized im. 
virtually any study involving human behavior: in. compensatory 


tracking tasks. 


В. COMPENSATORY TASKS AND QUASI-LINEARIZATION 

The compensatory tracking task, shown. im Figure 1,, 
assumes that the error signal is the only information that: 
the operator is receiving. In this study the operator 
attempts to minimize a visual error signal. by using: a hand- 
operated controller. Tracking situations such as this. are. 
often encountered in aircraft flight control; е.д.., a pilot. 
attempting to maintain some desired pitch attitude in. the 
presence ot — a taria terres... IC has Deen. dhovrw that. 
in tasks such as this, the operator is nonlinear and. time 
variant in behavior. He may, however, be successfully 
modeled in a quasi-linear fashion [Ref. I]. This quasi- 
linearization implies that his response to visual stimula- 
tion is largely linear and time invariant; i.e., his 
dynamics are largely those of constant-coefficient linear 
differential equations. To account for nonlinear and/or 
time-varying behavior, the model includes a remnant signal 
as shown in Figure 2. The remnant is that portion. af the 
operator's output which is not linearly correlated. with the 
input. The human operator model thus consists of a linear 


describing function, ох determined from the quasi-linear 


10 





analysis, and a remnant, n(t). 

It should be noted that this quasi-linear: model is: of: 
little use if the remnant is relatively large, since them. 
the operator's behavior is predominantly nonlinear and/or- 
time varying. 

To determine the human operator model it is necessary 
to calculate the linear describing function Y, (34) ,. and the 
remnant power spectral density, Ó ,)س(‎ from physical. 
measurements oi Signals of finite duration in. a laboratory. 
experiment. The input signal must appear: to the operator” 
to be random, although it need not be truly: randon,, and. 
the operator must be well trained; i.e., nat undergoing. 
adaptation or learmung- ВеЁ. 21. 

In order to measure the human describing functiomn. and: 
remnant, one of three techníques can be employed.. These. 
are the direct Fourier analysis of the system signals;. the 
use of crosscorrelation methods, and a model optimization 
technique. These methods are discussed in Ref.. 3.. 

If, as done here, direct Fourier analysis or. cross— 
correlation methods are employed in measuring the human 
describing function and remnant, then the concept: of 
spectral analysis must be introduced. The next section is 
devoted to this topic. A more thorough treatment. can. be 


found in Ref. 4, Appendix D. 


Tdi 





II. SPECTRAL ANALYSIS 


A.  PERIODIC SIGNALS 
A periodic signal, f(t), with a fundamental frequency 
Wy and period T, satisfying the Dirichlet conditions. [Ref.. 4,, 


p. 579], may be represented by a Fourier series 


£(t) = 7 Е(пш)е? 15 


n=% 
where Ernie ~ | ae, 
-T/2 


The autocorrelation function for the above periodic 


signal is defined as 


ёт) = [ү (ttt) at. 


T Jaja 


о сап be weilitten 


шэг? тү "ыг 
£f Fu 


and is equivalent to 


$ ee 0) |} eena tà) f ecem ac E 
n=-© БЕ 


With F(n) denoting the complex conjugate of F(n), 
со 


фе (1) = V Е(0)7е)991" pu) - y р(п)р(п)е1П®1Т 


n=- œ т=—со 


ФЕ (т) = |0) Ь А 


n--o 
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The power spectral density, 9eg(n), is defined 


o m 1722: eein t ат 
f т ЕЕ 
-T/2 
and it can be shown that 
" | 2 
Ф--іп) - Е (0) Р 
It can be seen from this relationship that 


сз 


фрее(т) =) Фуга 


n=- 


jno, t 


The crosscorrelation function, фр ғ (Т), of two: periodic 
П 
signals may be found in a similiar manner if both. signals 


have equal fundamental frequencies, 0 and both. signals 


e 


satisfy the Dirichlet conditions. Assuming these. conditions 


are met, then 


, pil 
s +f £ (t)f (ttt) dt .. 
тыры т }_туз ї 2 


In a fashion similar to that for single signale, it can be 


shown 
6, Е, a 5 Р. (mE, (2) 0191 E н 
uu n--o 
The cross-power spectral density, E = (п), is defined 
2 
ЊЕ с = Td h. p, C) 5389 Ta, 
1 -T/2 


and it can be shown that 


ф (n) = F, (n)F., (n) 
ty sd % 2 


Using this relationship 


ф (т) = 90 т" | 
11, пса 12 
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В. TRANSIENT SIGNALS 
A signal, f(t}, is defined to be transient r£ 


im f(t) — O. 


ta 
If the transient signal, f(t), satisfies the Dirichlet 


conditions in amy finite interval, and if 


f 159 ШЕ deo , 


then the signal may be expressed as a Fourier integral [Ref.5, 


p. 279]. Under these conditions, the. Fourier: integral,. 


ЕЕ) = >| AE 


“--ОО 


gives the values of f(t) at all points, including those where the: 


fue tiIon 1s not cómtiñudus. The Pouriér transform of f(t) is: 


Т т--СО 


F (jo) -/ £ (ке Ese 


The autocorrelation function for the nonperiodic: signal 


is defined 


Фе (т) -fT Ee Etat x 


са 


This Can also be written 
ф--(1) - A " 
Letting | 
Ф-г-(0) = МЕТ |2 i 
where 9 ec (0) is defined to be the energy spectral density of 


tne Signal f(E), jt can be shown that 
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теат f jme ax 


00 
Thus it can be seen that the energy density spectrum and 
ehe autocorrelatáon function of a transient signal are az 


Fourier transform pair, 


| EM A = Jam 
and 


о е (ш) ADEL 


If two transient signals, f) (t) and f(t), each satisfy 


the Dirichlet conditions in all finite intervals,, and if 
y | £, (E) jdt < о and MIS " 


then 


- y _ BL IE = T x Jor 
te e 10) -/ Е, (Є) Е, (Є+т)ає = af Far, Gwe dus 


am OO 
Now with 


Ф (0) = E, (OYE) (qo) , 
EJES y 2 


where Š < Ё (ш) is defined as the cross-energy spectral. 


1 2 
density of the signals £, (t) and Е (Е), 1t can be shown that 


Ф (w) -f ф Бар ан | 
E -ө 5152 | 


and has as its inverse transform 


l ч JWT + 
ф E 2- | Ф (ш) е dw 
£1£5(1) = 27 k... £,£, 


The transient signals £, (t) and Е (Е) are said to be 


linearly uncorrelated when $e £ (т) =“ Олово at. 
` | 


b5 





С. RANDOM SIGNALS 
In general, a random signal, f(t), from a stationary,. 
ergodic random process [Ref. 4, p 279], does not: have: az 


Fourier transform since 
со 
Т | 566) | ёс 


is not finite. An autocorrelation function. may be. defined. 
for the random signal f(t) as 
З pun = £^ Ға (+T idi 
Toe мт 
Since ФЕЕ (Т) satisfies the Dirichlet conditions for: all 


finite intervals and 
Jemen s 


it can be represented by a Fourier integral. It car: then. be 


shown that 
^o NM ` шт 


where 9 eg (9) is defined as the power spectral density of 
the signal f(t). 9 ee (0) is the Fourier transform ОЁ фе (Т), 
e -јшт. 
$ ce (0) one dt 
If there exist two random signals, f. (t) апа £,(t), 
which are sample functions from two different random processes, 
each of which are stationary, ergodic, and. jointly ergodic, 
then the crosscorrelation functiom is given. as 
^ do 
1 | 


T 

| т Ї ха 

ф (т) = LIM 28 f. of (terrat = | ф (u) ^ ag, , 
£,£, Too 27 J-r 1 2 2т 1. 515 
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where $e f (ш) 15 defined as the cross-power spectral 
1 2 
density of the random Signals, [> (cy and £ (ts. 0. -- (m) has: 
: 2097 бв, 


the Fourier transform 


Ф (ш) -f ф zT - 
fif, 03 £,£, De CET 


The random signals £1 (t) and £^ (t) are. said. to. be 


linearly uncorrelated if ФЕ Е (=), = @ for all T.. 
1-3 
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III. DESCRIBING FUNCTION AND REMNANT RELATIONS 


A. FREQUENCY DOMAIN RELATIONS 
Again consider Figure 2 with the input a sample: function 
from an ergodic random process. It is seen. that 
Ее) = зш) - Со) , 
where the Fourier transforms areas defined. in. Section BE... Now 
c(jw) = [N(ju) + + GuJEGO] ү (бә) „ 
then 
1 (ju) - N(ju)Y, (du) 
5029) = ру Gay Gey | 
Fe] 
Finally, after multiplying by I(jw), the complex: conjugate 


of the input, | 
IGo)|rGo) - wGo)Y, Go)| 


I (joe GOY = | ; | 
|- + Eli 0009 





In a similar manner, 
P(jo) = TAMEN + М(јо) 
and then substituting for E(jw) from above, 


је бег (о) : N Go) Y. Go) Y. Gu) 


P(jw) = “арт хүрэл зэр ee 


апа B 
Lc | T (ju) Y, Go) ro) -- N Go) 
РР = AAA 
р E 
Likewise, T(jw)C(ju), E(ju)P (ju), E(jo)E (jo), E(jw)P (ju), 
and € (Ju)C(Jw) may be calculated. The results, to be 


utilized shortly, are shown in Table I. 
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В. FINITE RUN LENGTH 

Since it is impossible to have experimental. runs of: 
infinite duration, measurements using finite time histories: 
are necessary. Reference 6 indicates: that: finite: run. 
lengths can be handled analytically as follows:: If i (t) 
15 the input and defined 


LL G) 


inet) = | 


HS w ST 


f 


0: ELSEWHERE. 


7! 


then this function can be considered to be- transient- and- 


have a Fourier transform 
I (jw) -| i (ie аа 


The other system signals and their transforms can be defined: 
in precisely the same manner. If the run time, I. DNE 
enough to ensure accurate power spectral. measurements,. yet; 
finite so that the respective Fourier transforms exist, 


then the following spectral relations are valid [Ref..6]. 


®..(w) = 21Ҹ hi [roo]? 


ir To TAM 
and | | 
_ LIM I -,. x 
$a Mu == Lk Towel 


where 9. , (w) and $, (o) are power and cross-power- spectral 


> 


dems1tres respectively. Utilizing Table r, 
1 1 (30) 1 (30) Y, (10) + І(јш) М(јш)Ү (39) | 
ыг MI XN AAA је) ај 


LIM 


T>% 


Pip (0) = 





p' 
but 


I EP | 
LIM | ı  Тбојт(јој|ј= Ф.. (о) 
го Ба 29) 2 ¿34 10) 


189) 





and 


ат Ез u o) |= ЭС) 
thus 
Y (jura; y (w) +O, (w) Y (Ju) 
Pu Е i ша ок. 
Now, 
è (w) = Ј mia К 


Since by definition the remnant, n(t), is linearly. uncorrelated: 
with the input, i(t), then PO = 0 for all т.. Thuss 
2.2, (9) - 0 апа 
Y Mw) Ф. .:() 
Фр ШТ = р Чан | 
хул чэн нэ 
In a like manner, 


== І(1ш)Е(3ш) x г’ 


DEMO. LL m! 
X cbe E 
эг 
у ҮЭ Ф: (ш) > 0, (ш)Ү (Ju) 
le 


1 + Y Go)Y. (ju) 


Again since n(t) and i(t) are linearly uncorrelated, 


mU) = 0 for all t, then $, (o) = 60. Thus: 


Ф (0) 
Фей. = er _ er 37 
те 


PAIS 


Also | 
"LIM ІШ ари С 
Тоо ۶رمن ۴ ا‎ tao) 


aed i 18! BE (430) EG) rom (jo)+N (jo) ]) 
pp T>% |2T e. | | 
| 1+ СЗАО l + о | 


um 4 
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Or 


Е 2 | Ex 
w $1; (0) | Y4 Go) | +@ (0) Y, (0) +9, (0) Y, (3) +9, (6) 
РР 


ЧЕ? 2 
Ley (ju)Y..(ju) 


Again, since Ф to) $ ; (o) = 0, 


Ñ 2 
- 222407 x Go) | +8, (0) 


PP |Ley, Gu) Y, Gu) | 


C. DESCRIBING FUNCTION AND REMNANT 
From "ЭОД 2.2 0), апа vae Y uim) and 0, (0) may 

be found. -Utilizing Ó, (س)‎ and Y шш 

2-04) = 916 (0) [L+Y, (Gu) Y, (o) 
and 

21-23 Фф... М AY o (TAD Y (Tu) 
ғаш т) їр | р С | 
р 

Thus, 

. " Í | iuga HER. Ж 

Фу. (9) l+y (jo)Y (je) = a бар s S o. sa E 


ЈЕ 


Or 
- P on 


Qu UO) 


У 410) 
E Jet 


Also, from Ф (ш), 
РР 
ó (um) — |14Y (jo)Y. (ju)F^8.. (u)-1Y. (30) | ^9, (a) 
nn p c ‚ РР | p | li 
In addition, Y. (Ju) can be determined and calculated. from 
Pre (ә) = $n (0) Y. (ju) + Bigs gu Y. Cua) 


Again $. (w) = 0, thus 


zl 





2.- (5) = AS E шш) T 
But 


o (4) т, (і) = cm) p 


thus 
л) 
22 (ја) = NON .. 


The functions in Table IT form the: basis: of: thesdeseribing 


function measurement technique: utilized in this: study. . 


р, SINUSOIDAL INPUTS 

The relations in Table II. are predicated on.the. existence” 
of a random input. In experimental work, а random appearing: 
input is often used and can be generated as a summation of 
sine waves [Ref. 3]. The input can. thus be represented by: 

n 
i =% A, sina, t 
k=1 | 
where the u, аге chosen such that in. a. finite. run length. 
there will exist an integral number of periods: or: complete 
cycles. In this analysis a run time of 150 seconds was: 
used and 0.08 < ú, < 40.0 radians per second. 

Utilization of a sinusoidal input results in. system 
signals that have both random and periodic components. With 
what is now a mixed signal, the question. arises as to which 
power spectral relationship should be. used.. The. solution. is 
to use the periodic power spectral relationships: for-measure-- 
ments made at the input frequencies and to use the. random, 


finite relationships at all other frequencies. 
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It should be noted that if the experimental run: length, 
T, is large and contains an integral number: of- periods: of: 
each of the input sinusoids, then the Fourier transforms: of: 
the periodic and finite random signals differ: only’ by. az 
constant of proportionality [ВеЕ.. 6]... 

In the periodic case, it was shown. that. 

цаа 1 11 Ex. TR. 
F, (r) = f Е Е т. 


2T 
kJ Ur 


and for the random signal, 


F (Jw) 4 Fr (t) MET 


If it is recalled that Esp (E) = (0 for t<@f" and teT,, amd. 





letting T = M Thr where m, is the number: of periods, Тр, ОЁ 
Егеачепсу Qu; then 
Ж š i ІШЕ Б... 2 М 
F, (n) = z- f К (је ЈАЧЕ = = f E KE ea ta, 
m, T 
k -Тү КК ~m, Ty 


or 
T Е -10, Ё 
Е. (п) = AJ f(t)e J"x^at  .. 
k 2Т ! 
-Т 
thus it Can be seen" that 
2TF, (n) = E (Ju, ) 


It should be noted that the same expansion. of: the: limits 


en the integral cam be used for the periodic: function;; thus. 


Y ff x: —jnu, t- 
Prem) = arf, uf serene ктат , 
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or 


Ян ша ue is 
o--(n,) = SH f(t)f(t+r)e k ‘atar 


Alsea Гос tre fim e candem Function. 


ti 
Ж | je Á 


Or 


= | WT rar s 
ФЕ (Ш) = x. f, F(t) £, (ttt) e ;; 


then by equating the integrals, since £, (t) = ite for 


НТ, Be TSNSSen tat at a specific input frequency,w,, 


(w = 2T 


Ф ғ (О) т gg Pp) r 


where 
а me. T : 2: 
Pee ye = مب‎ Fa | | ES 


ЇЕ has been shown that 


шог 


“ae (oy 


Y (ju) 
Р о: 


At the input frequencies, with T large and containing an 


integral number of periods of each input frequency, 


Y (ju) = "grec a арт = En I 
p É 2 9 ie dy r E 
Or " 
er a FS _ Ey) 
k 9; (nk) Ein )Ein) Пк) 
Similarly, 
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а ей š С (пу) 
€ k AU Р (пу) 

This illustrates that the cross-power spectral measure- 
ments need not be made and only the Fourier transforms are 
needed. The latter is usually an easier measurement than the 
former. 

The terms in 2а (9) can be examined in a similar manner. 
In this case, if measurements are taken at frequencies, Wn r 
other than those used in the input, then the expression for 
Ф да (9) is somewhat simplified; i.e., 

ап (Әһ) = |1*X, ТМ ИМ) 
since Ф.. (0) = 0. 
It should be emphasized that at frequencies, Qu, other 


than those used in the input, 


| abor) | Од 
ӘЗІРШЕ > ЕТІНЕ)! апа Y. (39) Pin) , 


since I(n, ) = І(п,) = 0. Thus BU, and Yo (Jup? must be 
estimated, as direct calculations can be made only at the 
input frequencies. In order to estimate ШЕГІН simple 

linear interpolation between Хэг and рт) сап Бе 


used. 
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ТУ. COMPUTER MECHANIZATION 


А. CEXPERTMENTAL SET=UP 

The measurement of a human's describing function and 
remnant in the compensatory tracking task of Figure 2 was 
made using a hybrid (analog-digital) computer. The error 
signal, e(t), was viewed as the vertical displacement of a 
horizontal line on an oscilloscope screen. 'The operator's 
controller consisted of a non-moving force stick. Control 
was effected by fore and aft pressure on the stick; e.g., 
if the line on the oscilloscope moved above the datum, the 
operator applied forward pressure to move the line down, and 
vice-versa. The input, i(t), and controlled element dynamics, 
у. (Јо), were mechanized on the computer as were the measure- 
ment algorithms to be described. Each experimental tracking 


run lasted 150 seconds. 


B. FAST FOURIER TRANSFORM 
The pertinent relationships are again 


| P (ny) 


and 
Pan (09) = [14% Guy Y, Go) ۶۵ (4p) | 
Now 
Р (пр) = A- f * s ie etae ; 
-T/2 
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МЕ 


Л T/2 -Г/2 
P (n,) E y if p(t) cos (nu, ©) dt + if pit) sin(nu, +) GE 
-T/2 -T/2. 


These integrals may be approximated. by the following: 


summations: 


. АЕ Ы E а А а 
P (ny) = —= 2 P(nAt)cos (w, nat) Бар ^p 2, p(m qe sin (a, nat) 
ТЕ 
N 
Pr =2p (nAt) cos (w, nAt) 
n=0 
and 
B N 
Pr =), p(nAt) sin.(w.ndt) „ 
n=0 i 


then the fast Fourier transform P (n, ) can. be. written: as. 
ams ЖЕ ЕК к: 
РОТЕ s %, 3| 


Similarly, it can be shown that 





cur TA. FUE E mU DNE ES e 
C (nj) = — | сы F and E. (n) — ЖЕ | e. o, | 
From this, RE 
| . EX Pk 
у Јер) д + JB 
ёр e 
s k 
and 
E 2 | 
а, 9 12 
Е (ар) | = – AA 
Р дов. 
NTC К: 
It can also be shown that 
В. B 
бий - Uk — Px 
TN = 2--22 | Бүү Е = 
$ о Јер) Байт R- tar” 
К Pk: 


2” 





in order to validate the simulated controlled element- 


dynamics, у. (је), on-line measurement of Y. (ju) cam. bes 


utilized: A 2 E Zu 1/2 
=. E "C 
Үс (Зе) | = nz 
! A +B 
Ex Ex. 
and 
Вр Во. 
а. E. -1 12 -1 
$ Y Go) — tan. ne tan. M 
Ру: Сү. 


The power spectra of the remnant can also be: determined. 
from the above relationships with the interpolation. process: 
described earlier. The determination. of the power: spectra 
of the operator output, n ОКУ can be accomplished: using 
the measurements p(t) at any desired frequency [Ref..90].. 


Previously it was shown that 





8 ww Y ó e = Gran) o ; 
thus 
ò (w) > P(n )P(n )2T = 2T|P (n. ) | ШИШ, DES ET 
pp h 107” wu | h. | e Ч Ph E 
and finally, 
o fu y о зен ја, | 


From this it can be seen that in order to determine. the 
Operator's describo tunctidwr, (im J” the remnant, 


Pan Mp) and the controlled element dynamics, Y, (Ju the 


к), 
only measurements needed are those of the error; operator: 
output, and controlled element output. If each of: those 


measurements, taken at specific times, is then multiplied 
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by the proper trigonometric function and. summed over: the: 
entire run, then the describing functions and remnant car. 
be calculated. 

The major drawback of the fast Fourier transform 
technique is the necessity of calling the trigonometric: 
functions sin(x) and cos(x) during the run. This requires: 
so much computer time that the multiplication. and addition 
computations cannot be performed between analog-to-digital 
interrupts. This in turn means that all data. must: be stored 
for later computational purposes. 

The necessity of calling trigonometric functions during 
the run can be avoided if the relationships for: sin(a+b) and. 
cos (a+b) are utilized in recursive fashion. If the initial 
time of the run is "a" and the time between measurements is 


"b" then 


sin (a+b) sin (a)cos(b) + cos(a) sin. (b) 
and 
cos (atb) = cosí(a)cos(b) 4- sin(a)sin (b). 

It is obvious that a recursive process can be mechanized 
which obviates calling trigonometric functions during the 
run. 

The use of this method yields results immediately upon 
completion of the run and conserves computer: storage: space. 


Within seconds of run completion, the data may be analyzed 


from either numerical or graphical read-out.. 
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W. RESULTES. 


The operation of the program was checked by measuring 
the transfer function of know elements or filters in. 
place of the operator. The results are shown. in. Figure. 3;,, 
ц, 5, amd 6. 

In validating the remnant measurement: technique,, the. 
operator was again replaced by an element with Known. transfer: 
function. A signal 


SW 


n(t} = Asin (wt), Wy 


io e" 
was inserted at the output of this element.. Measured 
values of NC, were then compared with. the theoretical 


A? 
— value [Ref. 7]. 


Ч 

Actual describing function and remnant measurements. 
were then taken on two subjects. The first: af these: had 
considerable previous experience. The second subject had 
experience only with this experimental set-up.. The measure- 
ment results are shown in Figures 7-18 for controlled 
elements of Ү (5) = 1.0 апа =. Each of these figures 
represents rhe results of EU tracking runs. 1 а ишта 
the human describing function contains the gain of the 
controller. 

The describing functions illustrated are comparable 
with those obtained by other experimenters;: e.g..,, [Ref.. 8].. 

This computer program represents a powerful tool. for 


use in experimental investigations involving a human 


30 





controller. It can, for example, be utilized imn a variety of” 
situations in which quantitative models of pilot. behavior- 

are desired. The program itself requires little core storage.. 
This means that considerable storage is available for: simu— 
lating complex aircraft dynamics, providing detailed display 


formats and calculating performance measures.. 
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| + Y (jo) Y (j o) | + Yo) e Yd jo 
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— 


Ng ш) + Elj Y Gu Ie Gu] ING) + El o) Yg Gu) [Ye w ] 


32 





TABLE IL 


Ys C o) ne a 





Y (| ш) = ic 
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mue) = Фе) ПУ о Дал = (al 


DW) Yo? + Ф (о) 
Фр = | F: Y (je) Y (jo ) |° 
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